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I. INTRODUCTION 

In this paper we consider a question similar to that 
posed in the title of Ref. [T]: How much does a sphere 
rotate when rolHng on a curved surface? In Ref. [T], the 
old problem of the rotation of a torque free, non-spherical 
body is reanalyzed. The angle of rotation is identified to 
have two components, one dynamical and one geometri- 
cal (the so called Berry phase), independent of the time 
elapsed during the rotation. Here we consider a related 
but different problem: a sphere is made to roll without 
slipping on a given curve F on a surface. The question 
is, if the sphere completes a circuit, what is the rotation 
matrix connecting the initial and final configuration of 
the sphere? The problem we are considering is therefore 
a kinematic rather than a dynamic one: the trajectory 
of the contact point of the sphere and the surface is dic- 
tated externally and the rolling constraint is imposed. 
We make contact with recent approaches that consider 
the same problem [51 [3] (but on a plane), in particular, 
we address a nice question posed by Brockett and Dai 
[4]: a sphere lies on a table and is made to rotate by a 
flat plane on top of it, parallel to the table. The question 
is: if every point of the plane describes a circle, what is 
the trajectory and motion of the sphere? 

We treat the problem by exploiting its isomorphism to 
the precession of a spin 1/2 in a time-dependent magnetic 
field. In the mapping, the arc length of the curve plays 
the role of time. For rolling on a plane the magnitude 
of the magnetic field is 1/R with R the radius of the 
sphere, and the direction of the magnetic field is that of 
the instantaneous angular velocity of the rolling sphere. 
For a curved surface the normal curvature and the torsion 
of the curve affect the value of the effective magnetic field. 
Closely related to the present paper is the use of of the 
isomorphism between classical dynamics and that of a 
spin 1 /2 by Berry and Robbins in Ref. [5] , especially their 
classical view of the Landau-Zener "7^ problem. From 
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a pedagogical perspective, the novel contribution of this 
paper is to use the isomorphism to discuss rolling spheres 
on an arbitrary surface. 

The precession of a spin 1/2 is widely treated in the 
literature and one can borrow those results to acquire 
an intuition for the rolling sphere. Conversely, since a 
rolling sphere is a tangible physical problem, the present 
treatment can be useful pedagogically in presenting spin 
precession, Berry's phases and it's classical counterpart, 
Hannay's angle [B]. 



II. ROLLING ON A PLANE AND QUANTUM 
PRECESSION 

Consider a sphere of radius R rolling on a curve F 
on a plane. We define a local triad of unit vectors at 
the contact point (the so called Darboux frame JF): the 
tangent t to F, the normal n to the surface, and u = nxt, 
the tangent normal. For rolling on a plane n is a constant 
vector, and the velocity of the center of the sphere is along 
the tangent to the curve. This situation will change for 
rolling on a curved surface, but, as we will see, the general 
idea of the mapping to a precessing spin is the same. 

The translational velocity of the sphere is V = iV{t) 
and the rolling constraint means that the instantaneous 
velocity at the contact point is zero [9]: 

TS X (nR) = V = tV{t) (1) 

with uj the angular velocity and R the radius of the 
sphere. This equation is nonintegrable and constitutes a 
paradigmatic nonholonomic constraint [10 . 

Taking the cross product with n on both sides of the 
above equation we have 

V(t) V(t) 
To = — X t = — ^u. (2) 
R R ^ ' 

Notice that in the above equation we have used the "no 
spin" condition • n = 0, that is, we are consider roUing 
without an instantaneous rotation along the normal. 
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The instantaneous velocity X of a point of coordinate 
X (with respect to the center of the sphere) on the surface 
of the sphere is 



x = Tzr xx = 



m 

R 



u X X. 



(3) 



Now we rewrite V{t) — ds/dt where s is the arc length 
of the curve T{t), and ([3| becomes 



dX 

ds 



u 

R 



X X. 



(4) 



If we regard X — (x, y, z) as a magnetic moment, the 
above equation describes its precession in the presence of 
a magnetic field B — —ji{ux,Uy,Uz) = —To of constant 
magnitude 1/R. The direction of B is — u, and varies 
varies with s, the arc length, which plays the role of time. 
If the rolling is on a horizontal plane, then B^—O, but we 
keep this notation to make contact with the rolling on an 
arbitrary surface. 

There is an isomorphism between the rolling sphere 
written in this way with a spin 1 /2 precessing in this 
magnetic field. This isomorphism can be seen clearly if, 
(using B = —~uj) we rewrite Equation (MJ) in the form 



(5) 



which is the same as the following equation of motion for 
two complex numbers a and b (we write s instead of t for 
time in order to keep the analogy) 





ds 



B 

Br + iB. 



y 



Bx — iBy 
-B, 



(6) 



with the identification 



X = ah* + ba* 
y = i (ab* — ba*) 
z = aa* — bb* . 



(7) 



The real numbers (x, y, z) represent the coordinates 
of a point on the surface of the sphere referred to a 
coordinate system fixed in space (that is, not rotat- 
ing), and whose origin is in the center of the sphere. 
The above mapping is certainly possible because of the 
SU{2) - SO{i) isomorphism [U,. 

Equation ([6| is Schrodinger's equation for the spinor 
X = {a, b) in the presence of a magnetic field B: 



Equation ([t]) implies that we can extract the behavior 
of the rolling sphere as a function of arc length by solving 
the motion of a spin 1/2 in a time- varying magnetic field. 
To our knowledge the equivalence between the motion of 
rigid body and a two- level system (a spin 1/2), in the 
form of the mapping of Eq. ^ was first pointed out by 
Feynman, Vernon and Hellwarth [12 and later discussed 
several times [15]. Earlier, Bloch [13j had derived the 
precession equation for the density matrix of spin 1/2 and 
therefore the points {x, y, z) that result from the mapping 
from spinors are called the Bloch sphere. 

The pedagogical novelty of the present paper (an al- 
ternative title of which could well have been "The rolling 
of the Bloch sphere" ) is to discuss the rolling using the 
arc length as time and identifying the isomorphism be- 
tween the rolling sphere and the quantum spin in exactly 
solvable cases. 



III. WARMUP: CONSTANT MAGNETIC FIELD 

Consider the simplest case of constant magnetic field. 
We choose B = -Bok, constant in the +z direction. This 
corresponds to the sphere rolling on a vertical plane. Eq 
^ becomes: 



ds 



with solutions: 







'Bo 










. 


-Bo j 





a(s) 
b{s) 



(9) 



(10) 



Replacing ( 10 ) in Q we obtain 
'Bos 



x{s) = a;(0)cos 

y{s) = y(0)cos 
zis) = z(0). 
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Bos 
2 



y(0) sin 
a;(0) sin 



Bos 

2 
Bos 

2 



(11) 



which means that the sphere is rotating clockwise around 
a constant axis in the z direction. This corresponds to a? 
in the —z direction. In other words, a constant magnetic 
field in the z direction corresponds to the sphere moving 
in a straight line in the xy plane, rolling on a vertical 
wall. The same situation applies if a constant field is 
directed in any other orientation. 



IV. THE LOLLIPOP AND THE PLANAR FIELD 



as 



(8) 



where h — 1 and H the Hamiltonian. Also, the vector 
S = |(o'a;, CTy, CTz) is the spin operator, and ct^ are Pauli's 
matrices. Notice that in this mapping, the magnetic 
fields and the frequencies have units of inverse length. 



Consider a magnetic field varying on the xy plane as 
B = i? (cos as, sin as, 0). This corresponds to u rotat- 
ing with the same frequency in the same plane, and 
the rolling problem becomes that of a sphere of radius 
R = l/B rolling counterclockwise on a circle of radius 
r = 1/a (see Figure [T]). 




inxt = B_,,(i) 



FIG. 1: The lollipop, or a sphere rolling counterclockwise on 
a circle of radius r corresponds to a spin 1/2 precessing on a 
magnetic field that rotates in the xy plane. 



In turn, this corresponds to a time (or arc length) de- 
pendent Hamiltonian H = — B • S, which can be solved 
by noting that 



B • S ^ ( ^ )-U*(y^]U, (12) 



ith 



U 



giQs/2 g 

e-^""/2 



(13) 



Substituting the above relations in ^ we obtain a time 
independent equation for the coefficients xi-^) — {'^^ ^) — 
(e'"'/2a,e-*"''/26) 



ds 



a B 
B -a 



= H 



(14) 



Transformations (12 1 and (13 1 correspond to trans 



forming to a frame that rotates with angular velocity 
a [18]. When transforming to the rotating frame, the 
angular velocity acquires a component a — \/r in the z 
direction and the frequency of rotation in the rotating 
frame is 



rn 



i?2 



(15) 



This can be se en i n the spinor language by noting that, 
since H in Eq. (14) is time-independent , the solutions 



a B 

B ~a 



[cos{ns/2) 



m 

ia ■msm{ns/2)]x{0), (16) 



with ±ri = V-B^ + ck^ the eigenvalues of H and m a unit 
vector in the direction a/B = R/r. Equation (16 1 de- 



scribes a rotation at a rate f2 with respect to an axis in 
the direction of the "stick" of the lollipop (the direction 
joining A to the center of the sphere (see Fig. ([T])). No- 
tice that solving for the evolution by exponentiating H is 
possible because H does not depend on s. If there is an 
s-dependence and the matrices H at different s do not 



commute the solution is a "time ordered" exponential 
that in general is not exactly solvable. 

After the lollipop completes a circle, the angle 5 of 
rotation is 



2n, 



-n = 27r 



(17) 



Notice that, when R ^ r the angle of rotation is (5 = 
27rr/i?, corresponding to rolling in a line of length equal 
to the perimeter of the circle. 

In anticipation of the next section we mention than in 
this case, since the rolling is on the plane, there is no 
geometric phase. When the rolling is on a curved surface 
the situation changes. Notice that we are using the term 
"geometric phase" in its relation with the spin problem 
in the adiabatic approximation. This phase is different 
from the nonholonomy when the sphere of arbitrary ra- 
dius describes a loop. 

We see that, after traveling on a circle the sphere is 
rotated by 27rf2/a with respect to an axis tilted with 
respect to the plane; this is the nonholonomy treated in 
[3J and [17,. 

When the sphere rolls on a plane, and on a circle of 
radius r much larger than its radius R, it comes back 
rotated around an axis that lies on the plane, by an angle 
given only by the dynamical phase. The extra term that 
originates in the curvature of the surface is what we call 
the geometric phase. 

The angle of rotation S (of both the spin and the lol- 
lipop) has a simple geometric interpretation: when the 
lollipop rolls, the point of contact C moves on the cir- 
cular rim of the cone ABC (see Figure [T]) . At the same 
time, the point C "paints" on the sphere a circle of di- 
ameter BC = 2rR/^/r'^ + i?^. (This is easily calculated 
with simple geometrical considerations from Figure [l]) 
This means that after a revolution of length 27rr the an- 
gle rotated is 2nr / {BC /2) from which Eq. (17 1 follows 
immediately. 

At this point we consider Brockett's question men- 
tioned in the Introduction. Notice first that, as the 
sphere rolls on a circle, the velocity at the top of the 
sphere is twice the velocity V at the center of the sphere. 
Since each point of the plane on top of the sphere de- 
scribes a circle of radius i?i , the velocity Vp of the plane 
also describes a circle. Therefore, since the sphere has a 
rolling condition with the upper plane, then Vp = 2V, 
meaning that, as the plane describes a circle of radius Ri 
the sphere describes a circle of radius Ri/2. 

We showed this with a nice classroom demo: on a piece 
of paper draw a circle of radius 5 inches (twice that of 
a tennis ball). Orient the label of the tennis ball at 45 
degrees with the vertical (the sphere is going to roll on a 
circle of radius r = R, and therefore the axis of rotation 
is going to be at 45 degrees and the precession frequency 
will be, from (nsl, V2). Paint a mark on a transparent 
glass, which in turn will serve as the upper plane. Also 
mark three points on the circle separated by /3 = 127 
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degrees (7r/-\/2). Looking through the glass, guide the 
mark on the glass over the circle on the paper, and notice 
that, each time the glass rotates by (5, the tennis ball 
rotates by tt with respect to a moving axis at 45 degrees. 

Notice also that for s = 27r/a the spinor x changes sign 
due to the 1/2 factor in the transformation. Nevertheless, 
since the mapping of ([t]) is quadratic in a and 6, changing 
their signs corresponds to the same values {x,y,z) for 
the orientations. More specifically, the quantities a and 
b determine univocally x, y and z, but the reverse is not 
valid: the quantum evolution determines univocally the 
classical evolution but there is some ambiguity in going 
from the classical to the quantum case. For example if we 
perform the "gauge transformation" (a, &) e'^'^^^\a^b) 
the mapping to the X coordinate remains unchanged. 

We will come back to this point in the next sections 
when we discuss the geometric phase for rolling. 



V. ROLLING ON A CURVED SURFACE 

In this section we extend the treatment of rolling on 
a plane to rolling on a curved surface (See Figure |2]). 
If we call Xp the coordinate of the contact point, the 
coordinate Xj, of the center of the sphere is: 



Xc = Xp + i?n, 
and its velocity is given by 

(in\ ds 
R- 



(18) 



ds I dt 



(19) 



The rolling condition is that the velocity of a point of 
the sphere in contact with the surface is zero (See Eq.(|l])): 




FIG. 2: Sphere rolling along a curve F of zero torsion (mean- 
ing that the velocity of the center of the sphere is parallel to 
the tangent of the curve at the contact point). 



The discussion for the planar case extends to the 
curved surface, and the rolling of the sphere is equiva- 
lent to a spin 1/2 precessing on a magnetic field B(s) 
given by 



B(s) 



(1 — K„i?) U -t- Trt 



(24) 



with the arc length s playing the role of time. In the 
following section, as an example of this formulation we 
consider rolling on a spherical surface. 



VI. SPHERE ROLLING ON A SPHERICAL 
SURFACE 

In this section we consider a sphere of radius R rolling 
on a second sphere of radius r. The rolling line will be 
a parallel of latitude it /2 — 9 (see Figure |3|. This means 
that the normal curvature is constant 1/r, and also that 
the torsion is zero. The magnetic field for the correspond- 
ing spin problem is therefore: 



~uj X (ni?) — Xc. 



(20) 



Again, taking the cross product with n on both sides 
of the equation above we obtain 



— n X Xc 
R 



(21) 



We now replace ( 19 1 in (21 1, and use the fact that, for 



a curved surface, the variation of the normal is given by 

-K„t — T^U, (22) 



ds 



with K„ the normal curvature and the torsion of the 
curve, both evaluated at the contact point. We obtain 



(1 — KnR) U + Trt 



ds 
'dt' 



(23) 



B(s) = - 



(25) 



with R± = rR/{rzLR) a reduced radius and the plus and 
minus signs refer to the rolling outside and inside of the 
sphere of radius r respectively. 

For a sphere rolling on a parallel, the instantaneous 
angular velocity (and the magnetic field) describes a cone 
forming an angle 6 with the vertical. The total arc length 
of the parallel is r sin 9 meaning that the vector u rotates 
with angular frequency a given by a = l/(rsin0). The 
corresponding magnetic field is therefore 

B(s) = {Bx,By, Bz) — ^i— (cos 6* cos as, cos 6* sin as, — sin( 
R± 

(26) 

with the term B • S in the corresponding Hamiltonian 
given in this case by 



B S = 



1 1 



COS ye 



sin 9 cos 9e 

icts 



sm( 



(27) 
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Notice that, if we compare with the rotation in a plane 




FIG. 3: Sphere rolling on a sphere. 



This again is an exactly solvable Haniiltonian that was 
first studied by Rabi 



Using the same transformation matrix of Eq. ( 13 1 the 



above Hamiltonian can be rendered time independent. 
We write it in the following form 



from Eq. ( 15 ), the "outer" rotation corresponds to rolling 



on a circle of radius equal to that of the unfolded cone 
tangent to the parallel (See Fig. [s]). The angle of rotation 
along that circle is not 2n but 2ticos0. This geometric 
factor is the same that appears in Foucault's pendulum 
and in Berry's phase for a spin precessing on a cone (we 
will come back to this point below). Also, notice that 
when r = R the angle of rotation is always 2n indepen- 
dent of latitude. 

The "inner" roll has the interesting feature that, when 
r = 3i? the angle of rotation is 27r regardless of latitude. 
This amusing feature can be verified easily at the equator: 
roll a penny inside of a circle of radius three times the 
radius of the penny and verify that the penny completes 
a full rotation in the rotating frame (and of course two 
full rotations in the lab frame) . 

We finish this section with a discussion of the differ- 
ences and similarities between the Berry phase for a pre- 
cessing spin 1/2 in the adiabatic approximation and the 
rolling of two spheres. 

The Hamiltonian for a spin in a magnetic field that 
precesses along the z axis at frequency a is given by ( 28 1 , 



where in principle a and Bq are independent parameters. 
If a ^ i?o (the adiabatic approximation) the eigenvalues 
(eigenfrequencies) of H are 



fv _ 1 / -Bosin6' + Q! Socos6' , , , 
2^ Socos6' Bosm9-al' 

with Bo = l /R± ■ 

The eigenvalues of H are 



2aBQ sin 6* ~ Bq 



(33) 



After a period of time 2n/a the change A0 in the phase 
of the spin is 



2Ra 



R 



± 



rsmO 



(29) 



with the spinor precessing, in the rotating frame, around 
an axis that forms an angle (3 (see Figure |3| with the xy 
plane, with 



tan (3 = tan 6 



R 



1 



r + R sin 9 cos ( 



(30) 



The second term in ( 30 1 reflects the fact that the small 



sphere rotates instantaneously on the tangent plane that 
contains BC (see Figure |3]) . Equation ( 30 ) can be easily 
derived by simple geometric considerations from Figure 

After a complete revolution the angle or rotation 5 is 



5± = 2'Krsm6U,±. 
After a little algebra we obtain 



(31) 




27rWl-|-sin2 6l— ' 



R^ 



(32) 



A(j) = 2tt 



Bo 



2'!rsm0. 



(34) 



The first therm is the dynamical phase and the second 
is a purely geometrical one, independent of the parame- 
ters Bq and a, and give by (half) the solid angle described 
by the field. 

For the rolling sphere we can also study an "adiabatic 
approximation" since a <^ Bq corresponds to r ^ i?. 
In other words, in general the adiabatic approximation 
will correspond to the radius of the rolling sphere much 
smaller than the radius of curvature of the surface. On 
the other hand, in contrast with the spin case, the fre- 
quency of rotation a = 1/r sin "knows" about the 
latitude and the curvature. So we expect some differ- 
ences and some similarities. Re placing the values of 
Bq = 1/ R± = 1 /R ± 1/r in ( 34 1 we obtain the angle 



of rotation of the sphere in each case (in the adiabatic 
approximation) 



A(j)± ^ 27rrsin6' ( — ± - I - 27rsin( 
R r 



r sin 9 

27r^^ 



(Outer rolling), 
47r sin 6 (Inner rolling) 



(35) 
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The above interplay of curvatures for inner and outer 
roUing is a special case of more general treatments of 
kinematics of rolling and is discussed in Ref . [52] . 

From Eq. (351, we see that in the outer rolling case 
there is no Berry phase, something we could have ex- 
pected because of the analogy with the rolling on a flat 
plane. The angle of rotation is in this case given simply 
by the rotation on a straight line of length equal to the 
perimeter of the parallel. However, for the inner rolling 
we indeed have a geometric phase twice as big as that 
of the spin 1/2. Our treatment is a nice example of the 
appearance of a geometric phase in a classical system, 
originally discussed by Hannay 0. 

In the next section we discuss the general connection 
between rolling and the Berry phase for spins in the adi- 
abatic approximation. 



VII. THE ADIABATIC APPROXIMATION AND 
ROLLING ON A CURVED SURFACE 

In this section we compare the equivalence between 
the adiabatic approximation for a spin precessing in a 
magnetic field that changes direction at a slow rate and 
rolling on a surface. In the spin case, the dimensionless 
parameter controlling the approximation is the ratio of 
the instantaneous frequency (proportional to the instan- 
taneous magnitude of the field) with the rate at which 
it's direction is changing. 

In the rolling case the instantaneous frequency corre- 
sponds to the magnitude of B(s) and the rate of change 
in its direction is related to the normal curvature and to 
the curve's torsion. 

In the adiabatic approximation for spins |20j . one 
works in an "instantaneous" basis, treating first s (time) 
as a parameter and solving the eigenvalue equation as 
though the problem were static: 



H{s)x{s) = n{s)x{s). 



(36) 



Then the general solution is written as linear combina- 
tions of the instantaneous eigenstates. As a result, in the 
adiabatic approximation, the spinor at time s is given by 



(37) 



The argument of the second exponential above repre- 
sents the dynamic phase, which involves the integral of 
the following angular frequency: 



n{s) = \B{s)\ = y[l-n^{s)R]' + [r{s)R]^ 



1 

R 



K„(s) 



(38) 



This can be seen, for example from Equation (27): the 



eigenvalues of B • S with s treated as a parameter are 
±|B(.)|. 



The (instantaneous) direction of the field is in the di- 
rection ub given by 



B(s) (1 -K„i?)u+r,t 

— n-,^M — , (,o9J 



ms)\ 



(l-K„i?)" + r2 



In general, the eigenvalues of a Pauli matrix in an ar- 
bitrary direction ub • o given by the unit vector = 
{ux,Uy,Uz) are ±1. This is verified by noting that (defin- 
ing Ux + iuy = pe"^) 



(ub • <t) X±(ub) 



X±(Ui3) = ±x±(ub), 
(40) 

with x±(us) = (1,±(1 - u^)e'^"^/p). Notice that the 
dependence of x on s is through the orientation of u. 

The first term, the geometric phase 7, is the Berry 
phase, and is given by 



7(s) 



ds'x{uB{s'))^^x{nB{s'))- (41) 



If the rolling describes a complete circuit, 7 measures 
the solid angle described by u^. This can be seen ex- 
plicitly as follows. If we express in polar coordinates 
Ub = (sin 9 cos 0, sin 9 cos (p, cos 9) then the spinor in that 
direction is: 



X(us(s)) 



cos 



(42) 



This means that X^^X = the inte- 

gral over a closed circuit F can be written as 



7r 



A • dl 



(43) 



with A = §^u^. Since V x A = — u^, using Stokes 
theorem, the line integral of A is the flux of a monopole 
in the origin, giving the solid angle |23j . 

Notice that this solid is traced not by the normal to 
the surface but by u^. This means that the solid angle 
measures a combination of the normal curvature and the 
torsion of the curve. In contrast, the solid angle traced 
by the normal measures the geodesic curvature |24j . In 
summary, we have shown that, when a sphere of radius R 
rolls on a surface of local radius of curvature and inverse 
torsion much larger than R, the angle of rotation S in a, 
closed curve of length L is given by 



dsKn{s) — S, 



(44) 



with S the solid angle traced by u^. Note that, if we 
specify this result to the sphere rolling on the parallel of a 
sphere of radius r, we have L = 27rrsin0, k„ = ±l/r (for 
inner and outer rolling respectively) and S = 2t: sin 9. 
Replacing these in Eq. ( 44 1 we obtain the result of Eq. 
(35 1 as expected. 
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